Two-Particle Correlations in the Wave Function and Covariant 

Current Approaches 

Dmitry Anchishkin^'Q and Ulrich Heinz^'P 

^ Bogolyuhov Institute for Theoretical Physics, 03068 Kiev, Ukraine 
^ Physics Department, The Ohio State University, Columbus, OH 43210 
O ■ (Dated: February 2, 2008) 

o ■ 

(N 



(N 



X 



We consider two-particle correlations, which appear in relativistic nuclear collisions 
due to the quantum statistics of identical particles, in the frame of two formalisms: 



O 

Q 

ly^ ' wave- function and current. The first one is based on solution of the Cauchy problem 



whereas the second one is a so-called current parametrization of the source of sec- 
ondary particles. We argue that these two parameterizations of the source coincide 



^1^' when the wave function at freeze-out times is put in a specific correspondence with 

j-^. a current. Then, the single-particle Wigner density evaluated in both approaches 

^ [ gives the same result. 
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I. INTRODUCTION 



Q , The models and approaches which are used to describe the processes occurred in the 



reaction region in relativistic heavy-ion collisions are examined by comparison of provided 
predictions with experimental data on single-, two- and many-particle momentum spectra, 
which contain information about the source at the early stage (photons, dileptons) and 
at the stage of so called "/reeze-ouf (hadron spectra). Two-particle correlations or the 



Hanbury-Brown-Twiss interferometry (HBT) encapsu 



ate information about the space-time 



structure and dynamics of the emitting source Usually, consideration 

of the correlations, which occur in relativistic heavy-ion collisions, assumes that: (i) the 
particles are emitted independently (or the source is completely chaotic), and (ii) finite 
multiplicity corrections can be neglected. Both approximations are expected to be good 
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for high energy nuclear colhsions with large multiplicities. Then, correlations reflect a) 
the effects from symmetrization (antisymmetrization) of the amplitude to detect identical 
particles with certain momenta, and b) the effects which are generated by the final state 
interactions of the detected particles between themselves and with the source. On the first 
sight one can regard the final state interactions (FSI) as a contamination of "pure" particle 
correlations. But, it should be noted that the FSI depend on the structure of the emitting 

n 

source and thus provide as well information about source dynamics [8|. 

Several surprising questions motivated by new experimental data appeared recently in the 



HBT. For instance, the experimental measurements on two-pion correlations jol. [lol. Illl. Il2| 
give the ratio of Ront/Rside ~ 1, what is much smaller than that predicted theoretically (the 
so called "RHIC HBT Puzzle"). This raises the question to what extent some of the model 



predictions are consistent with experimental measurements |13l . |14| or may be the observed 
discrepancies are due to such an "apples-with-oranges" comparison. All this drew attention 
and inspired a more detailed discussion of the theoretical background of the HBT. In the 
present paper we are going along this line, we would like to clarify a question concerning 
different kinds of parametrization exploited in the HBT. 

The nominal quantity expressing the correlation function in terms of experimental dis- 
tributions [2! is 

where Pi (k) = Ed?N/d?k and P2{KM) = EaEb d^N / {d^kad^h) are single- and two- 
particle cross-sections. 

In the absence of the final state interactions the theoretical expression for the two-particle 
correlator reads 

2 



C{q,K) = l± 



J d^X e^^-^Si^X, K) 



d^XS{X,K I d^Y s(y,K-^ 



(2) 



where K = (fc^ + kb)/2 , q = ka — k^. This expression was obtained in the different 
approaches. In the so called "wave-function" approach source function S'(X, K) is defined 
in the following way 

S-wf (X, K) = j d^x e^^-^ Yl P77' ^7 + f ) i^Y - f ) ' (3) 

7>7' 
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where p^-^i is the density matrix which in thermal equihbrium has the form p^-^i = 
S^Y exp {—Ey/T). The wave function, ip^{t,x.), is taken at freeze-out times, i.e. t G ts- 
Freeze-out hyper-surface E is a spatial surface which moves in space in the same way as, 
for instance, the surface of the balloon during pumping. It represents an imaginary border 
between two domains: inside the surface a strong dynamics takes place whereas outside 
the surface the particles propagate outward freely. Wave function at freeze-out times can 
be regarded as initial one for its further history and because its further evolution is free 
(we do not discuss final state interactions so far) it can be easily taken into account. As it 
intuitively understood the free evolution can be reverse back and resulting cross-section and 
other measurable physical quantities, for instance source function S{X,K), are determined 
through initial values of the wave function, i.e. by the values of the wave function at freeze- 
out times. Rigorous evaluations give exactly this result. On the other hand, the strong 
dynamics which acts inside freeze-out hyper-surface results in creation of the quantum state 
at freeze-out times. Hence, the wave function at freeze-out times is a final state of the strong 
dynamics. Representing experimentally measured quantities with the help of these states 
we can study strong interactions in dense and hot nuclear matter. Because of this creativity 
the separation of the interaction scales in space and time which is made with the help of 
freeze-out hyper-surface looks so attractive. 

Correlation function (121) was derived first in the model where essential point is a 

n n 

parametrization of the source by use of the currents J-y{x) [ij (see also [I5j) which become 
then the constituent elements of the source function 

5cur (X, K) = J d^x e*^-^ Yl Pii' Ji{x + l) j;, (X - f ) . (4) 

7.7' 

As a matter of fact, both approaches should give the same result in the region where 
they are valid. The goal of this paper is to find relation between source functions ([3]) and 
(jl]) obtained in wave function approach and covariant current approach respectively. 



II. SINGLE- AND TWO-PARTICLE CROSS SECTIONS WITHOUT FSI 

In this section we consider the two-particle quantum statistical correlations when one 
neglects the final state interactions of the detected particles. This phenomenon is visualized 
most transparently on the bases of the standard quantum mechanics. First, we briefly con- 
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sider the so called wave function parametrization of the source in nonrelativistic approach. 
This approach allows one to include also into consideration the final state interactions j^. 
Relativistic picture is considered on the base of the current parametrization and then on 
the base of the wave function parametrization of the source. First, we compare these two 
approaches in a non-relativistic sector and put in correspondence the source functions ([3]) 
and (jlj). After that the same comparison is carried out for relativistic sector. 



A. Wave function parametrization of the source. Nonrelativistic approach 

The probability to register two-particles which are created in the relativistic heavy ion 
collisions and have definite asymptotic momenta and k;, is compared usually with the 
probability to register independently two particles with the same momenta. That is why, 
we first turn to consideration of the single-particle spectrum. 

Let us consider a single-particle state emitted by the source. Its propagation to the 
detector is governed by the Schrodinger equation 

^^^^ = Mx)^,(x,t), (5) 

where /i(x) = — 2m^^' '^^^ index 7 denotes a complete set of 1-particle quantum numbers. 
Equation ([5]) is solved by ^/'^(x, t, to) = exp [—ih{ii.){t — to)] '?/'7(x, to) in terms of the single- 
particle wave function at some initial time to? see FigiH For the spherically symmetric 
fireball the values of the wave function '?/'^(x, to) parameterize the "freeze-out distribution" 
of the particles inside the sphere of the radius Ri as it is depicted in FiglH We assume that 
the detector measures asymptotic momentum eigenstates, i.e. that it acts by projecting the 
emitted single-particle state onto 0k'^*(x, t) = exp [zk ■ x — 'iti;(k)t] , where a;(k) = k^/2m. 
The measured single-particle momentum amplitude is then 

74^(k, to) = lim / ci^x0^"*'*(x, t) ■?/'^(x, t, to) . (6) 

t >OQ J 

The single-particle probability to detect the particle with certain momentum is obtained by 
averaging ([6]) and its complex conjugate with the density matrix p^y defining the source 
dynamics. This density matrix is characterized by a probability distribution for the single- 
particle quantum numbers 7 and by a distribution of emission times to. We write 

Pi(k) = P77' (k, t,) a; (k, to) . (7) 

7,7' 
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Figure 1: Sudden freeze-out hyper-surface FqFiRq for spherically symmetric fireball. 



We define the single particle Wigner density of the source which with accounting for emission 
times, t = to, reads 

S (^^^, k)=J d\x, - X2) - to) 6{xl - to) e^^'-(^-^^) ^ ^,{x,) ^;,(x2) . (8) 

7.7' 

This function accumulates all information about the source which emits the particles. Mak- 
ing transformation to new coordinates, X = {xi + X2)/2, the source function 
gets the form 

S{X, K) = 5{Xo - to) I d^x 5{xo) e^^'^ Pii' ^7 + f ) " f ) • (9) 



7.7' 



Then, the expression for the single-particle spectrum ([71) can be rewritten with making use 
of the source function 

Pi(K) = j d^XS{X,K). (10) 

Note, the factor 6{Xo — to)S{xo) in (Q carries information about space- like hypersurface 
where initial values of the wave function are given. For the sake of simplicity of the general 
scheme we start our consideration from a flat hypersurface, t = to, depicted in FigUl We turn 
to an arbitrary hypersurface in the next sections where relativistic approach is elaborated. 

Let us consider a two-particle state emitted by the source. Its propagation to the 
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detector is governed by the Schrodinger equation 

.9V'^(Xa,Xb,i(:) 



dt 



if(Xa,Xfc) V'7(Xa,Xfe,t) , (11) 



where -ff(xa, Xf,) = /i(xa) + h{-Xb). The index 7 denotes a complete set of 2-particle quantum 
numbers. Equation ffTTl) is solved by 

V^^(x„,Xfe,t) = exp[-iH{xa,yib)it - to)]ip^{yia,yib,to) , (12) 

in terms of the two-particle wave function at some initial time to- Detector acts by projecting 
the emitted two-particle state onto 0p"*p^(xa, x;,, t) = e*(P'''''""'^°*) e^^P'"'"''"'^''*) , where Ua,b = 
p^j,/2m^. We will only consider the case of pairs of identical particles, = mi, = m. The 
measured two-particle momentum amplitude is then 

A^(ka, kfe) = lim / d^Xa d^Xb 0k"X(^«' ^) V^7(^«' t) . (13) 

We assume the two particles are emitted independently, implying that at some freeze-out 
time ta the two-particle wave function '?/'^(X(j, x;,, t) factorizes 

V'7(x„,Xf„t„) = [i)^^{yia,ta)i)^S^b,ta) ± V^^,, (Xfc, t„) V^^Jx„, t„)] . (14) 

The indices 7a,7fe on the single-particle wave functions now label complete sets of single- 
particle quantum numbers. The time moment ta = to is the emission time of the latest 
emitted particle. Because of the symmetry of the wave function (IT^ it does not matter 
what time is nominated as latest one, ta or tb. By this we assume that symmetrization 
occurs when the last of the two particles is frozen out from a strongly interacting bulk. 

After hermitian inversion of the evolution operator and applying it to symmetrised (an- 
tisymmetrised) out-state two-particle amplitude (|T3|) gets the form 

where Xa = to and x° = to and by relabeling the variables of integration we transferred 
symmetrization from the state (IT^ onto out-state. By this we represent the measured two- 
particle momentum amplitude as projection of non-symmetrized two-particle wave function 
taking at emission times onto symmetrised (antisymmetrised) plane waves taking as well at 
emission times. 



7 



The two-particle probability to detect two particles with momenta and is obtained 
by averaging two-particle amplitude f|T5l) and its complex conjugate with the density matrix 
defining the source. This density matrix is characterized by a probability distribution for 
the two-particle quantum numbers (7a,7fe), also we average by a distribution of emission 
times {ta, h). We write 

P2{ka, kb) = J2 P7a7.' A„'76' (^a, kfe; to) {K, ^b, to) • (16) 

7a76:7a'76' 

We made the ansatz p^yaibaany ~ Piaia' Pibiy "which factorizes initial density matrix p-yaibaany 
in such a way that independent emission of the two particles is ensured. 

After straightforward algebra we write expression for the two-particle probability 

P2(g, K)= j d'XS (x, ^ + I) / d'YS (y, - I) ± 

± J d^Xe^'-^SiX.K) j d^Ye-"'-^S{Y,K) . (17) 

Finally we get the two-particle correlator C{q,K) ([2]), as a ratio of two-particle probability 
( IT71) and single-particle probabilities (ITOi) . where the source function S{X,K) is defined in 
accordance with Eq. ([9]), i.e. all integrations are taken at emission times or on the freeze-out 
hyper-surface. 

B. Current parametrization of the source 

Let us consider a single-particle state emitted by the source which we parametrized by 
the "source current" J^(x). Its propagation to the detector is governed by the Klein-Gordon 
equation 

{d^d^ + m^)^^{x) = J-,{x) , (18) 

where d^j^d^ = d"^ — V^. The index 7 denotes a complete set of 1-particle quantum numbers. 
(In a basis of wave packets these could contain the centers X of the wave packets of the 
particles at their freeze-out times t.) 

Single-particle momentum amplitude is defined as projection of the wave function "^yix) 
at "detector time t = x°" on the out-state 0^"*, 

^(k) = lim UT'*ix), ^-yix)) = lim / d^xcl)°i^''*{x,t)idt ^7(x,t), (19) 
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where by definition a{t) dt b{t) = adtb — {dta)b. We assume that the detector measures 
asymptotic momentum eigenstates, i.e. that it acts by projecting the emitted single-particle 
state onto 

t) = e-^-W^+^i^- = , (20) 

where a;(k) = V + . Then, momentum amplitude can be rewritten as 

A^{k) = hm / d'x f dStt^'\x) I Gr{x - y) J,{y) . (21) 

Substituting to ^ the Green's function {f'tr\x) = e'^^^'^^''^-'') 

Gn{x -y)=^e{x'- / ^^^^^^ [ft\x) ft^^\y) - f[r\x) ft^-{y)_ , 

which satisfies equation {d^d^ + m'^)Gii{x — y) = S'^{x — y), and using orthogonal proper- 
ties of the basic functions, / cPx fl"^^'* (x) i dxo fkt\x) = ±(27r)'^ 2co'(ka) 5^(ka — k;,), and 
/ d^x fj,'^^'* (x) i dxO fk^\x) = 0, after straightforward calculation we come to the answer 

A,{k)=zJd'ye^-^''y^"-'^-U,iy), (22) 

where integration is taken over infinite space-time volume and that is why the finiteness in 
space and time of the particle source which we deal with is accumulated in the "cut-function" 
Jj{y). Moreover, it should be pointed out that amplitude A7(k) in ( l22l) is nothing more 
as the on-shell Fourier transformation of the current, hence in this approach the amplitude 
to register the particle with certain momentum k directly reflects the model of the source 
which is settled by the particular definition of the current J^{x). 

The single-particle probability is obtained by averaging (j22l) and its complex conjugate 
with the density matrix defining the source. This density matrix is characterized by a 
probability distribution for the single-particle quantum numbers 7. We write 

i^i(p) = E PwMp)AyiP)^ (23) 

7:7' 

Inserting (1221) into (l23l) and using definition of the source function (jll) after simple algebra 
we come to result 

Pi(k) = J d^xS{x,k), (24) 

which gives the single-particle probability in the same form as we obtained in Eq. (ITOl) for 
the wave function approach but in contrast to the wave function approach the integration 
in is taken over infinite space-time volume. 
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Two-particle momentum amplitude is defined as projection of the symmetrized (untisym- 
metrized) two-particle wave function ilj^{xa,Xb) at "detector times x° — ^ oo and oo" 
where the index 7 denotes a complete set of two-particle quantum numbers, on the momen- 
tum eigen state 0k"*k6) 

A^{ka,'kh) = lim lim U^XAxa^Xb), '^j{xa,Xb 

= lim lim d^Xad^Xb(f)'^X,i^'^^^b) i dxO i d^,'} '^■y{xa,Xb) , (25) 

We label out-state by the values of measured momenta, i.e. and k;,. The out-state at 
detector times reads 

C!k.(^a,a;,) = /i:HxJ/i+)(a;,). (26) 

If the source is completely chaotic, i.e. the particles are emitted independently, that 
implies that the two-particle wave function is a product of two single-particle ones. For pairs 
of identical bosons (fermions) the two-particle wave function describing their propagation 
towards the detector must be symmetrized (anti-symmetrized). Taking the same arguments 
as for the wave function approach about delay of emission of one particle with respect to 
another one we write 

^lanA^a.Xb) = ^[^7„(Xa) ^7,(Xfe) ±^^^{Xb)'^j,{Xa)] = 
d^yad'^VbGR^Xa - Va) G R^Xb - yb)^[jia{ya) J-y^iVb) ± J-yaiVb) JytiVa) 



(27r)32a;(kJ (27r)32^(k,) ^"""^ ^""'^ M'^M Mh) ±J,Ah)Mka)\ , (27) 
where we use solution of the Klein-Gordon equation (fT8l) which is obtained with a help of 
the Green's function Gr{x — y) and 

lXk)= /Ae^'^('')^°-'^-%(2/), (28) 

is the on-shell Fourier transformed source current. We do not write in (!27|) the negative- 
frequency piece of the Green's function because it evidently disappears on the next step: a 
projection of the wave function \l/^„^^j^(xa, Xf,) onto out-state. Indeed, to obtain the momen- 
tum amplitude we substitute expression ( 1271) into ( l25i) and use orthogonality relations of the 
basic functions fl^\x). All this results in a simple final expression 
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The two-particle probability is obtained by averaging amplitude (129|) and its complex 
conjugate with the density matrix defining the source. This density matrix is characterized 
by a probability distribution for the two-particle quantum numbers (7a, 76) 

A,,,,(k„k,) A;,^^,(k,,k,) , (30) 

7a7b>7a'7i,' 

As in the wave function approach we made the ansatz P'^aibaany ~ Pi^u Pibibn which fac- 
torizes in such a way that independent emission of the two particles is ensured. 

Substituting momentum amplitude ( l29l) into ( l30l) . using definition (jlj) of the source func- 
tion S{Y,K)^ we can write for the two-particle probability 

P2{KM) = J d^XS{X, ka) J d^YS{Y,kb)± J d^Xe''^-^ S{X,K) J d^Y e-''^'^ S {Y, K) , 

(31) 

which coincide with expression (ITTI) obtained in wave function approach, consequently, we 
obtain correlator in the form ([2]). The integration on the right hand side of Eq. ( l3T|l is 
just taken over an infinite space-time interval, whereas in (ITTj) the integration is taken over 
freeze-out hyper surface, or over initial times. 



C. Wave function parametrization versus current parametrization. Nonrelativistic 

approach 

The goal of this subsection is to put in correspondence the "wave function" approach 
which was elaborated in the Section II. A to the "current" approach of the Section II. B. 
To do this we consider this correspondence first in the non-relativistic limit and then fully 
relativistic comparing will be carried out. 

We are going to obtain the current approach in the non-relativistic limit (see Appendix 
|A|) . We make a standard unitary transformation of the wave function to extract oscillations 
associated with particle mass \E'^(x) = e'"^"^^" ip-y{x) . With respect to the new wave function 
ip-y{x) the basic equation f|T8l) reads 

('^* ^ i " ' ^^^^ 

where 

j,(t,x) = --^e^'"*J,(t,x). (33) 
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and we skipped all terms of the order l/c? and higher, they serve as relativistic corrections. 
By this derivation we put in correspondence the current in the relativistic parametrization 
of the source with the current in the non-relativistic one. 

With respect to the quantum state ilj^{t,:s) the momentum amplitude can be rewritten 
in the following way 

A,(k) =Jim^/ d'xj d'yfi''^'*{x)Go{x-y)j,{y)=z J d^ e'^^'^'^y'-'^-^ 3,{y) = i~3,{k), (34) 
where 

is the Green's function which satisfies equation {idt + V^/2m) Go{x — y) = — y) and 
is the on-shell Fourier transformation of the current. 

Compare the amplitude flMl) . i.e. ^^(k) = ij^{k), with the correspondent amplitude 
obtained in the relativistic case ( !22l) we see that they coincide with one another, just in 
place of the capital letter J one should put a small one. Hence, the same transformation 
should be done in the definition of the source function (jl]). 

Non-relativistic Schrodinger equation (jdt — H{x.)^ "^-yix) = supplemented by the ini- 
tial condition, ip^{t = to,^) = $-),(x), can be written, as was shown in Appendix [B] (the 
generalized Cauchy problem [3]), in the form of the Schrodinger equation with the source 
on the r.h.s. of equation which is defined at initial time t = to, 

(tdt - i/(x)) ^^{t, x) = t <l>^(x) 5it - to) , (36) 

which is valid for t > to- Then, one can solve this equation with the help of the Green's 
function (135!) and write solution in the following form 



^p^{t-to,x) = i J d^yGoit-to,x-y)<^^{y) 



= e{t-to) I ^rfVe-*"(P)(*-*»)+^P-("-"'V7(to,x). (37) 

As a matter of fact, this solution coincides with that one obtained with the help of the 
evolution operator, '?/'^(t — to, x) = 9{t — to) e~*^*^P'''^^*"*°V7(^o, x), which we exploited in the 
wave function approach in paragraph III A[ 

On the other hand, solving the Cauchy problem in this way one can consider the expres- 
sion on the r.h.s. of eq. fl36l) as a specific current 

j^(t,x) = z7/;^(t,x)5(t-to). (38) 
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Then, going through all preceding steps to evaluate the source function S{X, K) with mak- 
ing use of this specific current one evidently discovers that the source functions (jl]) and 
([9]) coincide with one another. We regard this result as first example when two types of 
parametrization of the source can give the same answer for specific connection between 
current and wave function given at freeze-out times. 

We are going now to prove that the same is valid in more general case. Indeed, let us 
write once more the solution of eq. fl32|) 

?/'^(t,x) = j d'^zGoit - zo,yi- z) j^{zo,z) 

= i J d'^z d^y Go{t - to, x - y) ^0(^0 - ^0, y - z)jy{zo, z) , (39) 

where in the second line we split the Green's function at the point t = to using the group 
property of the Green's functions. We are making now the physical input: let us prepare 
the initial state, which will be used in the wave function parametrization of the source, 
^■yiy) = 4''y(to, y), in the following way (note, up to now we did not specialize a generation 
of the wave function at freeze-out times) 

$7(y) = J d^zGoito- zo,y-z)j^{zo,z). (40) 

Then, rewriting the second line in (l39l) with making use of the state ^^(y) (just defined in 
( HOj) ) we obtain the single-particle quantum state, ip^{x), at the times which are after t = to, 
i.e. after freeze-out, in the following form 

V^^(t,x)=i J d^yGoit-to,x-y)ip^{to,y) . (41) 

What is most interesting, expression (j^Tj) is exactly a solution of the Schrodinger equation 
(!36|) with i'-y(to,y) as the initial condition (see (1371) ). 

Let us make one note. In eq. fl39l) after splitting of the Green's function we meet the 
product of two ^-functions, 6{t — to)6{to — zo). Because we are interesting in detector times 
the value of time t goes to infinity and we can avoid the first ^-function. At the same time 
the second 6'-function cuts an action of the source current at the times t = to- But this 
feature does not distort the influence of the current if to is the freeze-out time. In other 
words, we assume that a life time of the current coincides with a life time of the fireball, 

j^(t,x) oc e{to - 1). 
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So, we obtain the same quantum state il)y{t,x.) in two approaches: 1) The current 
parametrization of the source, expression fl39l) . first hne, which is solution of eq. fl32l) : and 2) 
The wave function parametrization of the source, expression ( HTl) . which is solution of the 
Cauchy problem for the Schrodinger equation. If we start a description of the propagation 
of the particle to detector from the quantum state ^/'^(t, x) taken from (HTl) and go through 
all steps to the source function S{X, K) we come to expression ([9]) which is evaluated with 
the help of the initial states ?/'^(to,y)- In fact, as was shown in the section Hi At by this we 
obtain the source function exploiting the wave function parametrization of the source. On 
the other hand, we can start the evaluation from the same quantum state ip.y{t,x.) taking it 
in the form (|39l) (first line). Then, we come to the source function (jll) obtained in the cur- 
rent approach. Meanwhile, the starting point for both expressions is the same state '?/'^(t,x) 
(what gives the same amplitude, the same single-particle probability and so on). Hence, if 
the initial wave function, $^(x), and the current, j-yix), are connected to one another by 
eq.(l40l). then both evaluations of the source function S{X,K) give the same result. Thus, 
the single-particle Wigner functions constructed in both approaches are equal 

S{Y, K) = 6{Y' - to) J rf^y 5(y°) e^^'^ ^ {y + f ) {y - f) 

7:7' 

= / d'y e^^-y E P77' h {y + I) f,' {y-i)- (42) 

7.7' 

Consequently, the single-particle spectrum and two-particle correlations taken in the wave 
function parametrization of the source coincide with the respective spectra taken in the 
current parametrization. 

Connection fBOl) between current and wave function at freeze-out times has a transparent 
physical interpretation: the action of the current which describes in a semi-classical way a 
creation of secondary particles during the life time of the fireball can be accumulated in the 
wave function at freeze-out times. That is why, it does not matter what quantity is used then 
to describe the free propagation of the particles to detector. Moreover, the correspondence 
pij) . as it is seen in FiglH results in extension of the effective volume where initial wave 
function is given by adding a spherical layer for the radii r in the limits Rq < r < Ri. 
Indeed, all particles which were emitted during life time of the fireball from the boundary 
FiRf) accumulated now on the space-like segment F1F2. This means that the wave function 
'ip'^itQ, x) given on the extended space-like hypersurface, segment F0F1F2, takes into account 
all secondary particles which were "produced" by the source current. 
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D. Wave-function versus current parametrization of the source in relativistic 

approach 

First, we consider the wave-function parametrization of the source. Single-particle mo- 
mentum amplitude is evaluated as projection of the wave function, \E'^(t, x), taken at asymp- 
totic times onto out-state 0°"* 

A^{p) = lim / c?3x0°"*'*(t,x) i dt ^7(t,x) , (43) 

t >'CxD J 

where 0°'^*()f:, x) = e-«'^(p)*+«p-''. The wave function, \E'^(t,x) is solution of the Klein-Gordon 
equation, (9^9^ + m^)\E'^(x) = 0, which is supplemented by the initial conditions 

vl/,(t,x)|,^,^ = <|.,o(x), and ^Mli^ = $,i(x) . (44) 

Ot t=to 

As we show in Appendix [B] (see (1B17I1 ) this problem can be formulated as equation with a 
source which constructed with a use of the initial conditions (HH) 



{d^d'^ + m2)^^(t, x) = $^i(x) Sit - to) + <l'^o(x) 6'{t - to) , (45) 

which is valid for times t > to- Solving this equation with a help of the Green's function 
and inserting solution to (H3l) one can write the amplitude in the following form 



A,{p)=Jim^J d'x J d'y6{y'-to)f^-^^'*{x)t£^Gn{x-y)^%{y), (46) 

where /^+'(x) = 0p"*(t,x). Taking Gr{x — y) in the explicit form and using the orthogonal 
properties of the basic functions f^'^\x) we come to the answer 



^7(P) = / d'y5{y'-to) 



vl/W(to,p), (47) 



where '^{^\to,p) is the Fourier component of a positive-energy piece of the function 
^7(^0) y)- Note, the wave function consists from two contributions, positive- and negative- 
energy defined, ^E'7(a;) = "^^^x) + '^i^\x), respectively. Actually, for the sake of simplicity 
we consider a flat space-hke hypersurface, t = to, the segment F1F2F3, as it is depicted in 
FiglJl At the same time, the amplitude P7|) can be expressed in the covariant form as well 

A,{p) = J da>^{y) /W'*(y) vl;^(^) . (48) 
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Figure 2: A sketch of the freeze-out hyper-surface (sohd curve) for a spherically symmetric ex- 
pansion of the fireball. 

where = cT(y) is the space-like hypersurface on which the wave function and its derivative 
are given. 

Formula P7|) gives a parametrization of the probability amplitude by the initial values of 
the wave function and its derivative at freeze-out times. The next steps are the same as in 
the section Hi B[ The single-particle probability is obtained with a help of the density matrix: 
Piip) = E7,y P77' ^^(p) A'(p) = E^.y P77' ^^^^'*(^o, P)^y H^o, p)- If we now represent the 
functions \E'^^^(to,p) as the Fourier integral 

(^0, P) = / d'x e-^P-" ^(+) (to, x) = 1 d^x - to) e'^'^ {x) , (49) 

and insert it to -Pi(p) we come to the standard expression of the single-particle probability 
Pi(k) = / d'^X S{X, k), where we define the source function 

Six, K) = 5(X° - to) / f/'x(5(x°) e^^-^ Pii' ^i^^'* + f ) ^ " f ) ' (^0) 

7:7' 

So, after definition the source function (ISO]) , which is constructed with the use of the 
wave function and and its derivative at freeze-out times, we are ready to compare this 
parametrization of the source with the current one. 

First of all let us mention that both parametrizations evidently coincide when one takes 
the special form of the current, J^(t,x) = d'if^{t,:>i)/dt6{t — to) + \l'^(to,x) 6'{t — to), which 
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is expression on the r.h.s. of eq. (H5i) . Then, we come to eq. (fT8|l which is starting point in the 
current parametrization of the source. We show now that the same is vahd in more general 
case (our consideration is very close to that one developed in the section Hi CI) . 

The single-particle amplitude at asymptotic times (12T1) which we obtained in the section 
111 Bl can be written in the following way 



A^{k)= lim I d-^x I d^ye{x^ -y^^ 



(51) 



where we use representation flB20p of the retarded Green's function, Gr{x — y) = 
9 (x° — y^) G^~^\x — y) — G^^\x — y) , and orthogonality relation of the basic set of func- 
tions fk^\x) and fl \x). To distinguish the "current" amplitude ( l5Ti) from the "wave 
function" one (H71) we marked it by tilde. On the next step we use the group property of 
the functions G^+\x - y) (see ( 16221) ) 



G(^\x-z)^G'^+\z-y) 



Inserting this expression into (|5T|) and making the following convolution 



(52) 



d'x /i+)'*(x) t ^ GW(x -z)=^ ft^'*{z) , 



(53) 



we get the amplitude flSTl) in the form 

A^{k) = J d^yd^z 6(^z^ -to) 

d^yd^z 6 (^z° - to) 



d 



Jiiy) 



■My) ' 



(54) 



where the last line in flM|) is obtained under assumption that the source current "works" 
just during the life time of the fireball, i.e. J^(t,x) oc 6{to — t). Taking into account this 
feature one can define the wave function at freeze-out times as: 

^^{to, z)=J d^y Gnito - z - y) J^(/, y) . (55) 

Inserting this notation to the second line on the r.h.s. of (1541) one can rewrite the amplitude 
y4^(k) in the following way 



i^(k) = J d^z 6 [z° - to) 



_d_ 

dz^ 



(56) 
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Then, as will readily be observed the last expression coincide literally with the amplitude 
obtained in the wave function parametrization of the source fH7|) . Hence, we can write 

i^(k) = A^{k) . (57) 

Because, the group property can be written in covariant form as well G^~^\x — y) = 
/(^) d(T^{z) G^^^x — z) d /dz^"^ G^^\z — y), one can obtain the amplitude in the covariant 
form fj48l) . Then, equality fl57|) is valid for an arbitrary space-like hyper- surface. 

So, if we keep relation between current and wave function at freeze-out times in the form 
f l55|) we guarantee that the single-particle momentum amplitude will be the same in both 
approaches. The statement is valid also for two-particle momentum amplitude. Because 
the amplitude is the main constructive element of the single-particle probability ( l23i) and 
two-particle probability (15U]) . the equality of the amplitudes results in the equality of proba- 
bilities. This means that the source functions obtained in the wave function parametrization 
( 150|) and in the current parametrization of the source (jlj) are equal as well when eq. fl55|) is 
valid. 

It is necessary to clarify the time structure of the current J^(t, x). As a source of the 
single-particle state \I^^(t, x) the current acts during the life time of the fireball or when its 
time argument t is less than freeze-out times, t <ta: 

Jj{t,x) oc9{t^-\t\). (58) 

In the applications the cutting of the time interval is usually made in a soft way with a help 
of the Gaussian function, J^(t,x) oc exp (— t^/2r^), where r is of the same order as t„. Our 
previous consideration was based on the rapid cutting of the current on the freeze-out hyper- 
surface like that in (1581) . Can a smooth switching off destroy our scheme? It is necessary 
to point out that just the Fourier transformed quantities enter the single-particle and two- 
particle probabilities. Let us look at the shape of the Fourier transformed cutting profiles. 
It is interesting to note that the Fourier components of the both time cutting functions, the 
Gaussian function and 6'-function fl58l) . give approximately the same bell like dependence 
on energy variable E, J{E) = J^^dtJ{t) exp{iEt). These functions squared, J{E)^, are 
depicted in Figj3l Only a slight difference between these functions is seen and, therefore, the 
choice of the type of time cutting function does not affect much, at least qualitatively they 
give the same result. That is why, if we exploit the ^-function cutting rule (!58l) we obtain 
the same probabilities to register the particles as in the case of the Gaussian profile. 
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Figure 3: Comparison of the Fourier coefficients squared of the Gaussian and ^—function time 
cutting of a source current. 

III. DISCUSSION AND CONCLUSIONS 

We considered two types of a semi-classical parametrization of the source which give a 
transparent scheme of evaluation of the single-particle spectrum and two-particle correla- 
tions: the wave function and current parametrization of the source. The main ingredients 
of the wave function parametrization are the values of the wave function which are given 
on freeze-out hyper-surface (in relativistic approach the values of the wave function deriva- 
tive should be given as well). In describing a propagation of the particles to detector after 
freeze-out these values serve as the initial conditions in the Cauchy problem: all infor- 
mation about evolution of the fireball is accumulated in the single-particle wave function, 
■^7(^0; x), given at freeze-out times (see Figl2]). In relativistic case it is \l/-y(to,x) and its 
derivative d'^^ltQjX^/dt, then, the relativistic projection onto the out-state results that just 
the positive-energy defined part of the wave-function, ^^\tQ,x), is exploited. For the sake 
of simplicity we discuss here a flat space-like hyper-surface t = to =const. An arbitrary 
freeze-out hyper-surface is also considered in the paper. 

Once the wave function at freeze-out times, t = to, is given, then, the single-particle 
spectrum and two-particle correlations can be constructed with the help of the single-particle 
Wigner density S{X, K) which reads 

S{X, K) = 6{X^ - to) J d^x 6{x') e^^-^ Pii' ^1 + f ) ^7' " f ) ' (59) 

7:7' 

where the measure of integration appears as a result of the transformation: dl^xi d'^X2 6{x\ — 
to)5{xl - to) = d^X d^x 6{X° - to) S{x^) with X = {xi + X2) /2, x = (xi - X2) . To obtain the 
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source function in relativistic picture one should put in (159|1 the functions "^i^^x) in place 
of ip-yix), then, we come to expression (l50l) . 

We propose a scheme to generate the values of the wave function at freeze-out times. 
This can be done with a help of the current which parameterizes the source 

^^(to,x) = J dyU^yGn{to-y\^-y)J,{y^y), (60) 

where J^(?/°,y) oc 9{to — i.e. the life time of the current equals the life time of the 
fireball. If the wave function \l/^(to,x) and the current t/^(y°,y) which parameterizes the 
source are in relation ( !60|) . then, the single- and two-particle spectra evaluated in the wave 
function parametrization are equal to the same quantities in the current parametrization. 
This results in equality of the source function ( l59l) (or ( l50i) in relativistic approach) obtained 
in the wave function parametrization of the source with the source function (jl]) obtained in 
the current parametrization. Moreover, the correspondence (1601) results in extension of the 
space-like piece of the freeze-out hyper-surface by including the new piece which is created 
by the particles emitted from the fireball during its life, < t < to, through the time-like 
part of the freeze-out hyper- surface. As it is seen in FigH] the space-like part of the freeze- 
out hyper-surface, segment FqFi, is extended by adding a new piece, segment F1F2, which 
is created by the particles emitted from the boundary, segment FiRq. The same picture 
takes the place for the freeze-out hyper-surface of an arbitrary shape which is sketched in 
Figl2) the segment accumulates the particles emitted from the time-like part of the 

freeze-out hyper-surface during time span to. Hence, if the wave function at freeze-out 
times is generated by the source current as in ( |60l) . then, the wave function accumulates 
information about all particles emitted from the fireball. Moreover, the correspondence ( 160|1 
(the correspondence (j^Tl) in nonrelativistic it is seen in Figs JH [2] results in extension 

of the size of the fireball: the wave function parametrization reflects the radius of the system 
which is Ri (but not Rq). So, effectively the size of the system is bigger, an extended volume 
includes the spherical layer, Rq < r < Ri, which contains free particles emitted from the 
boundary of the fireball. 

However, there is a source of particles, for instance pions, which creates particles after 
freeze-out, for instance a decay of long lived resonances. It can be formalized by introducing 
of a "post freeze-out" current, Jpfo(x). Basically, from the very beginning we can separate 
current into two parts: the first part is a current before freeze-out, I-y{x), and the second 
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part is a current after freeze-out, /pfo(x), 

j^{t, x) = e{to - 1) i^{t, x) + e{t - to) ipUt, x). 

Then, the momentum amphtude (l56l) should be modified to the following form 

<-> 

A^{p)= J rf4x(5(x°-to)/i+^'*(x)z^vl>^(x)+^ |f/^x^(x°-to)4+^'*(x)/pfo(x), (61) 

where the second term on the r.h.s. of equation reflects the sources of particles which appear 
after freeze-out, i.e. for times t > tg- It turns out that these particles give contribution just 
into the single-particle spectrum. Because for the particles which are created after freeze-out 
a symmetrization, for instance of two-particle wave function, starts when the particles are 
separated by big distances the two-particle momentum probability, P2{K, q), has appreciable 
values just for small relative momenta, | q | < 10 -i- 20 MeV/c. These values of the relative 
momenta are not experimentally "visible" . 
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Appendix A: NONRELATIVISTIC APPROXIMATION OF THE 
KLEIN-GORDON EQUATION 

In this appendix we are going to obtain a relation between the current which appears as 
a right part (source) of the Klein-Gordon equation and a current (source) of the Schrodinger 
equation. For this purpose we take eq. ( |T8|) . which determines a current parametrization of 
the source, and consider it in a non-relativistic limit. This equation can be rewritten in the 
following way 

(9^,9^ + m^)^^{x) = -{idt + yjw? - V^) [idt - ^m? - V2)^^(t, x) = J^(t, x) , (Al) 

where t = x^. We make a standard unitary transformation of the wave function to extract 
oscillations associated with particle mass 

^^(x) = e-™^'V^(a;). (A2) 
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With respect to the new wave function ip.y{x) eq. (lAip reads 



{idt + m + \/m? — V^) (z(9t + m — V m? — \/'^)ip^{t, x) 



e^'"*J^(t,x). 



(A3) 



It is necessary to point out that from now on the energy operator idt with respect 
to the wave function ip^ix) is an operator of the kinetic energy because idtip^itjX.) = 
e*™* (i(9f — m) \[^-y(t, x). We shall also quote our consideration to positive enegies. That 
is why in the non-relativistic approximation we have the following inequalities for energy 
and momentum operators: {ihdt)/mc^ <^ 1 and {—ih'V)/mc <^ 1, where the broken brack- 
ets mean an averaging over some single-particle quantum state. Hence, the operator in the 
first bracket on the l.h.s. of eq. (\A3\i is a positive definite operator (it does not have zero 
eigenvalues). As a consequence, it always has inverse operator, that is why we write 



With making use of the relation — iV'?/'^(t, x) <^ m-?/'^(t, x) one can expand square root 
operators which we meet on the l.h.s. and on the r.h.s. of this equation in the Taylor series. 
Just keeping leading terms we arrive to the non-relativistic equation 



where we skipped all terms of the order and higher, they serve as relativistic corrections. 
A general scheme to obtain the non-relativistic equation of motion to any order of relativistic 
corrections from the relativistic equation in the presence of an external field was elaborated 
in Q. 

Appendix B: INITIAL CONDITIONS AS AN EXTERNAL CURRENT 

1. Cauchy problem for Scrodinger equation 

In this section we represent the initial conditions of a differential equation as an external 
current (the generalized Cauchy problem IigI). We consider a differential equation which 
contains a first derivative with respect to time. To be specific let us take the Schrodinger 
equation 



(jdt + 171 — ^vrfi — V^^ "ip-yit, x) = — (jdt + m + y/ — 



■) ' e*'"*J^(t,x). (A4) 




(A5) 




(Bl) 
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with initial condition: ilj{t = 0,x) = $(x). To solve eq. (IBip we define a new wave function 
x) by extension of ilj{t, x) for t < in the following way 



^(x, t) 



?/'(x,t), if t > 0, 
0, if t < . 



which obviously satisfies eq. ( IBll) if iIj{x) is a solution. As a next step we make the Fourier 
transformation of eq. ( IBll) separating the integral over time in two parts, dtF{t) = 
J'^'^ dtp it) + /o°° t;tF(t). Then, (EI]) reads 

rft e''^Hdti){t, x) - /7(x) / dt e^'^V(i, x) 

-CXD J ~00 

= ie*'^V(^,x) -i^uj dte*'^V(^,x) -i/(x)V(a^,x) = 0, (B2) 

where il){uj,yi) = J^^dte^'^^tjj{t,x.). After integration by parts one can insert the func- 
tion $(x), which represents the initial condition, to eq. (1B2I) keeping in mind that 
e*'^*^()f:, x) 0. Then, for the Fourier components we obtain equation 



(cu-if(x)) ^/'(cu,x) = z<l>(x). (B3) 
Let us make the inverse Fourier transformation of eq. ( 1B3[) with making use of the equality: 

— e-'-'uji,{uj,^)=idt / — e-^-*V^(u;,x) = 2 9,V^(t,x). 

-OO ZTX J-oo ZTT 

Then, one obtains equation in time-space representation 

[tdt - ^(x)) ^{t, x) = t <l>(x) 6{t) . (B4) 

The same equation is valid for the function ?/'(t,x). It is an interesting result because we 
reduce the initial condition for differential equation (IBip to impulse current which stands 
now on the r.h.s. of equation flB4p . 

One can define the Green's function (idt - H{x.)) G{x -y)= S\x - y). With a help of 
the Green's function one can write solution of eq. (]B4|) 

^Pit, x) = 2 1 A Git - to, X - y) ^(y) . (B5) 

On the other hand, it is solution of the Cauchy problem which was formulated as eq. (IBip 
with initial condition. Integral on the r.h.s. of flB5l) represents propagation of the initial 
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"excitation" $(y) which exists at time t = to to other spatial points x during time interval 
(t-to). 

In free case in the non-relativistic limit the Green's function Go{x — y) satisfies equation 
{idt + V^/2m) Go{x — y) = 5^{x — y), and it reads 



.mx^ 



(B6) 



where cc;(k) = k^/2m. Note, the non-relativistic Green's functions have the following group 
property 

j d'y Go{x - y) Go{y - z) = -i Go{x - z) . (B7) 

So, using the explicit expression of the free Green's function (]B6|) one can write solution of 
the Cauchy problem accumulated in eq. flB4p in the following form 



V^(t,x) = e{t) [ -f^ £ye-''^^''^'+'''<^-y^ ^y) = e{t) [ -f!^e-' '^■"$(k) 



.m(x - y) 



$(y) , (B8) 



2t 

where $(k) = / c/'^yexp (— ik ■ y) $(y). 

We turn now to another way of solution of the Cauchy problem flBlD . Formally one can 
write the solution of the problem as 

i){t, x) = e{t - to) e-^^(P''')(*-*«)$(x) , (B9) 

where p = —ihd/d^ is the momentum operator and 'il){t = to,x.) = $(y), in what follows 
we adopt to = 0. We are going to give the formal solution flB9p in coordinate representation 
and then compare it with (IBSp . Indeed, in coordinate representation (]B9p reads 

^(t,x) =0(t) /^^rfV(x|p)(p|e-^^(^^--)*|pO(p1xO$(x'), (BIO) 

where (x|p) = exp {ip ■ x) is the eigen function of the momentum operator. In free case 
the matrix elements of the Hamiltonian, i^o(P)X) = p^/2m, looks like, (p|e~*^"*^P'^^*|p') = 
(27r)'^5^(p — p') e~*'^''P^*, where uj{p) = p^/2m. Finally, for free case we can write flB9P in 
coordinate representation 

tPit, x) = e{t) [ e-*-(p)*+^P-(-— ')$(x') . (BU) 

We see that this expression coincides with (IBSp . So, we find that it does not matter in what 
approach one solves the Cauchy problem, with the help of the Green's function or using 
evolution operator, both approaches give the same result. 
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2. Cauchy problem for relativistic equation 

We consider the case of free scalar field which can be described by the Klein-Gordon 
equation 

(a^a^ + m2)^(x) = 0, (B12) 
which is supplemented by the initial conditions 

vl/(t = o,x) = $o(x), and ^^^i^^M = $^(x) . (B13) 

We are going to show how these initial conditions can be inserted into eq. (1B12I) as a specific 
current. As the first step let us extend the function x) (we define the function equals 
to zero for negative times) 



fr(x,t) 



^(x,t), if t > 0, 
0, if t < . 



We make the Fourier transformation of eq. flB12p separating the integral over time in two 
parts, J^^dtF{t) = J^^dtF{t) + J^dtF{t). Then, one obtains 

dt e'^^Xd^d^ + m^)^!{x) = / dt e'^'d^^{t, x) + (-9^ + m^) / dt e^'^*^(t, x) 

-OO J — CO J — CO 

^,^,avl/(t,x) OO r-^^^,^,5vl/(t,x) 



= -$i(x) + 2CJ<l>o(x) - LU^^iu, x) + (-(9^ + m^) ^{u, x), (B14) 

where ^(u;,x) = /f^ (it e*'^*^^(t, x). After two integrations by parts we insert the functions 
$o(x) and $i(x), which represent the initial conditions, to eq. (1B14P keeping in mind that 
e*'^*^(t,x) 0. Then, for the Fourier components we obtain the following equation 



(-cu^ ~dl + m^) ^(cj, x) = $i(x) - icj$o(x) . (B15) 



We make the inverse Fourier transformation of eq. flBlSP with making use of the equalities: 



r ^ e--(-..^) X) = d^t, X) , {-^ u) e"-* = ^ 6{t) . 

J-oo ZTT J-oo ZTT at 

Then, one obtains equation in the space-time representation 
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The same equation is valid for the function x) for times, t >0, 



{d^d'' + m^)^{x) = $i(x) 5it) + $o(x) 6'it) . 



(B17) 



So, we represent the initial conditions for the differential equation (]B12I) as a current 
which stands on the r.h.s. of equation ( ]B17[) . The use of the function 6'{t) is common: 



JdtF{t)6'{t) = -Jdt^6it). 



To solve eq.( ]B17l) we use the Green's functions Gr^x — y) (or Gf{x — y)). Because in 
further consideration we look for solutions \l/(a;'',x) at asymptotic times, oo, just a 

piece of the propagator Gr{x — y) which carries the positive defined frequencies really gives 
contribution. (Because of that it does not matter what kind of the Green's function should 
be used, retarded or causal one.) So, with the help of the propagator Gr{x — y) which is 
defined as, (9^9^ + rn?)Gii{x — y) = ^'^{x — y), we can write solution of eq. ( IB 171) 



^{x) = I d'yGnix - y) [$i(y) 5 + $o(y) 5' 
'dS5[y') 



- y) — — ^yy > y) 



(B18) 



The integration in flBlSP is going on the space-like hypersurface = 0. For arbitrary 
lypersurface a it can be written in the covariant form in the following way (see, for instance, 
2ol, ch. 7) 



^^^da^iy) Gn{x-y)^^<I^{y) 



(B19) 



with a^{y) as the space-like hypersurface on which the initial conditions "^{y) and d'^ly) / dy^ 
are given, i.e. these functions under the integral are defined when y E cr. By definition, 

hit) dt /2(t) = hmf2{t) - dtUl{t))f2{t). 

We derive now one more relation which is useful in the consideration. One can write the 
Green's functions in the following way 



Gr{x -y) = e(x''- [G^+\x -y)- G(-)(x - y) 



where 



G^^\x-y) = i J 



fi^\x)fr*{y) 



(±),*/ 



(B20) 



(B21) 



(27r)32cu(k) 

with fjf^\x) = e^^'^^ which obey the orthogonal relations. It is obvious that the functions 
G'^^\x — y) and G'^~\x — y) satisfy the Klein-Gordon equation, [d^d^ + m?')G'^^\x — y) = . 
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With taking into account normalization, / d?x f^^^'* {x)i dxO fp^\x) = (27r)^2u;(k) (5^(k— p), 
one obtains 



(B22) 

This can be written for arbitrary space-hke hypersurface a in the covariant form [2^ 



-y)= I da^{z) G^+\x - z) ^ ^^(z - y) 

J (a) OZ^ 



(B23) 
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